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In a heterostructure of graphene and the ferromagnetic insulator EuO, the Eu atoms induce 
proximity exchange and inter-valley interactions in the graphene layer. Constrained by the lattice 
symmetries, and guided by ab initio calculations, a model Hamiltonian is constructed that describes 
the low-energy bands. Band parameters such as proximity induced exchange splitting, spin orbit 
coupling, and inter-valley interaction are determined. Calculations of the Chern number identify 
the conditions under which the hetero-structures exhibit topologically non-trivial bands that give 
rise to the quantum anomalous Hall effect with a Hall conductivity of a^y ~ 2^ jh. 


I. INTRODUCTION 

Ever since the classification of the integer quantum 
Hall effect (IQHE) in terms of topological invariantsi, 
significant theoretical effort has gone towards realizing 
IQHE phenomenology at vanishing external magnetic 
fields. Haldane^ proposed that in the presence of an 
intrinsic spin-orbit coupling, spinless electrons hopping 
on a two-dimensional honeycomb lattice are topologi¬ 
cally non-trivial which can result in one-dimensional chi¬ 
ral modes along the edges of an insulator. Owing to 
the chiral nature of these edge modes, the gapless edge 
states are dissipationless and exhibit a Hall conductance 
(^xy = jh. This Hall conductance is a consequence 
of the Berry curvature associated with the Bloch bands 
in momentum space, and it is quantized only when the 
Eermi energy lies in the bulk band gap of the mate¬ 
rial. A number of other proposal have been made for 
realizing the quantum anomalous Hall (QAH) effect in 
mercury-based quantum wells^, optical lattices^, disor¬ 
der induced Anderson insulators^, magnetic topological 
insulators^ii, and ferromagnetic graphene^i^. The robust 
nature of charge transport which identifies the QAH ef¬ 
fect at vanishing magnetic fields might enable design of 
novel quantum devices for low-power electronics applica¬ 
tions. 

To realize the QAH state in realistic materials two con¬ 
ditions are necessary, (i) broken time reversal symmetry 
and (ii) topologically non-trivial bands. Since topological 
insulators (TIs) possess a large spin-orbit couplingiS, a 
route towards realizing the QAH effect is to introduce fer¬ 
romagnetic ordering in TIs. Immediately following recent 
successes in synthesizing magnetic TIs (MTIs), transport 
measurements in MTIs verified the predicted e^/h Hall 
conductanc o^d^ . Another approach is to engineer the 
QAH state in ferromagnetic graphene in the presence of 
Rashba spin-orbit coupling^. Recently, graphene was 
successfully deposited on an atomically thin-film insulat¬ 
ing ferrimagnet (YIG), and the transport measurements 
revealed an unquantized anomalous Hall effect due to 


proximity induced ferromagnetisnti^. Several other mag¬ 
netic material/van der Waals (vdW) materials combi¬ 
nations (for example graphene/EuO, graphene/BiEeOa 
and MoTez/EuOSdl-iS) have been proposed for possi¬ 
ble spintronic o^^d^ and valleytronicsi^ applications. In 
these systems, ferromagnetic ordering is induced by a 
proximity effect, which is an efficient way to make a non¬ 
magnetic two-dimensional crystal magnetic. Addition¬ 
ally, proximity induced magnetism by a magnetic insula¬ 
tor allows for control of the electron and hole densities by 
gating. Previous theoretical studies have reported an ex¬ 
change splitting gap of 36 meV in graphene/EuO hetero- 
structuresi^. 

In this paper, we investigate the possibility of a QAH 
effect in a EuO/graphene/EuO hetero-structure as shown 
in Eig. [T] The hetero-structure symmetries guide our 
construction of an effective Hamiltonian that captures 
the low-energy band dispersion. Ab initio calculations 
provide parameters such as exchange splitting and spin- 
orbit coupling. The model Hamiltonian captures the in¬ 
fluence of the Eu atoms which leads to two distinct types 
of inter-valley interactions in addition to the exchange in¬ 
teractions and Rashba spin-orbit coupling. An analysis 
of the graphene/EuO model Hamiltonian provides condi¬ 
tions under which the systems acquire a non-zero Chern 
number and exhibit a Hall conductance of 2e^/h. 

The rest of the paper is organized as follows. Sec¬ 
tion H presents the graphene/EuO hetero-structures, 
the density functional theory (DET) calculations, and 
the band dispersions resulting from the ab initio cal¬ 
culations. Section HI describes the model Hamilto¬ 
nian that captures the low-energy band dispersion of the 
graphene/EuO hetero-structures with and without spin- 
orbit coupling. In section IV, calculations of the Chern 
numbers identify the conditions which yield topologically 
non-trivial bands and the QAH effect. Section V con¬ 
cludes with a discussion of the possibility of proximity 
induced exchange and the observation of QAH effect in 
other graphene/ferromagnetic hetero-structures. 
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FIG. 1: (Color online) (a) Schematic view of a hetero-structure with graphene between two EuO layers. The O layers are 
terminated with H atoms, and the surfaces abutting the graphene are the Eu (111) planes, (b) The reciprocal lattice 
corresponding to the unit cells shown in (c) and (d) maps the K and K' points of the hexagonal graphene Brillouin Zone 
(BZ) indicated by the outer red hexagon to the F point of the commensurate BZ of the graphene/EuO unit cell indicated by 
the central green hexagon. The gi’s are the reciprocal lattice vectors of the hetero-structure unit cells in (c) and (d). 
Elevation and plan views of the unit cells corresponding to the two graphene/EuO geometries are shown in (c) for the 
Eu-misaligned structure and (d) for the Eu-aligned structure. 


II. FIRST PRINCIPLE CALCULATIONS 

The band dispersions of the EuO/graphene/EuO 
hetero-structures are calculated using the Vienna ab ini¬ 
tio simulation package (VASP)i^“— in the projected- 
augmented-wave methodic. The generalized gradient 
approximation (GGA) of the Perdew-Burke-Ernzerhof 
form^Si^— is used for the exchange correlation energy, 
and a Hubbard-U correction is used for the magnetic 
insulator, EuO. The on-site Goulomb repulsion and ex¬ 
change interactions on the Eu atom 4/ orbital are 8.3 
eV and 0.77 eV, respectively, and on the O atom 2p or¬ 
bital, they are 4.6 eV and 1.2 eV, respectively^^. The 
kinetic energy cutoff is 520 eV for all calculations. Dur¬ 
ing all structural relaxations, the convergence tolerance 
on the Hellmann-Feynman forces is less than 0.03 eV A. 
An 8x8x8 Monkhorst-Pack k-point mesh is used for 
bulk EuO. The calculated bulk lattice constant is 5.186 A 
which is very close to the previously published first prin¬ 
ciple calculations^^ and consistent with the experimental 
results. The lattice constant oq of graphene is 2.46 A. 
This results in a lattice mismatch of less than 1% in the 
3ao X 3ao unit cells shown in Figs. [IIc,d). 

The EuO/graphene/EuO structures consist of 
graphene between the (111) Eu planes of EuO. The 
hetero-structure with graphene on the Eu-terminated 
surface is more stable than graphene on the 0-terminated 



FIG. 2: (Color online) Band structure of the Eu-misaligned 
structure without SOC. Left inset: close-up of the 
low-energy band structure of the misaligned structure near 
r. Right inset: Band structure of the Eu-aligned structure. 


surfacei^. The relaxation of the 2D hetero-structures 
uses the same level of theory, cutoffs, and tolerances as 
described in the previous paragraph with a Monkhorst- 
Pack k-point grid of 4 x 4 x 1. A vacuum buffer space 
over 25 A is included to prevent interaction between 
adjacent slabs and hydrogen atoms passivate the outer 
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oxygen layers of the EuO films. The relaxed vertical 
spacing between the Eu and C layers is 2.517 A for the 
misaligned structure of Fig. me) and 2.555 A for the 
aligned structure of Fig. m- These distances are close 
the value of 2.57 A found previously for a single-sided 
heterostructure of graphene on EuO^^. 

Figs. IHcjd) show the two different hetero-structures 
considered. They differ by the alignment of the EuO- 
monolayer on opposite sides of the graphene layer. In 
both cases, graphene is placed on the (111) surface of 
EuO. This gives a commensurate hetero-structure with 
a lattice constant 2 times the lattice constant of a EuO 
unit cell and 3 times that of the graphene unit cell. In 
the aligned structure, shown in Fig. [T](d), the top EuO- 
monolayer is directly above the bottom EuO-monolayer, 
whereas in the misaligned structure shown in Fig[ljc), 
the top Eu-monolayer has an in-plane displacement of 
1.22 A with respect to the bottom EuO layer. In both 
structures the Eu atoms either sit at the center of the 
hexagonal graphene unit cell or at the the bridges of 
the C-C bonds coinciding with the inversion symmet¬ 
ric points of graphene’s honeycomb lattice. Therefore, 
in-plane inversion symmetry is preserved for both cases. 
However, as a result of the lateral displacement of the 
EuO layer in the misaligned hetero-structure, inversion 
symmetry perpendicular to the graphene sheet is bro¬ 
ken in contrast to the aligned hetero-structure where 
this symmetry is preserved. These symmetries play an 
important role in determining the band dispersion and 
the model Hamiltonian of the graphene/EuO hetero¬ 
structure. 


A. Band dispersion without SOC 

Fig. [2] shows the calculated band dispersion in the ab¬ 
sence of SOC for the misaligned structure of Fig. HKc), 
and the right inset shows the band dispersion for the 
aligned structure of Fig HKd). Both band dispersions 
are calculated along the path M — F — K of the com¬ 
mensurate BZ of the 37V x 37V graphene lattice. K = 
(27r/3a, 2-\/37r/3a, 0), M = and a = 7.38 A 

is the lattice constant of the hetero-structure unit cell. 
The most striking difference in the two band dispersions 
is the presence of a gap between the spin-resolved bands 
of the aligned hetero-structure, whereas in the misaligned 
hetero-structure the spin-up and spin-down bands in¬ 
tersect. The calculated values of gaps for both hetero¬ 
structures are tabulated in Table HI The energy band gap 
is denoted by Eq. The gap between the spin-up elec¬ 
tron band and spin-up hole band is A-j., and the gap 
between spin-down bands is Aj,. The spin-splitting of 
the electron and hole bands are 5e and 5h, respectively. 
In Table U and Fig. [51 the positive value of Eq = 127 
meV indicates a band gap between conduction and va¬ 
lence band, whereas the negative value of Eq = —38 
meV indicates a spin resolved band overlap. The values 
of A^ and A^ in the misaligned structure are half of their 


values in the aligned structure. Another striking feature 
of the calculated band dispersion is that the low-energy 
bands of the combined hetero-structures appear at F and 
have curvature. In contrast, to the ab initio studies of 
a graphene/BiFeOa hetero-structurei^, the Dirac cones 
are no longer at the K and K’ points, but at F, consis¬ 
tent with earlier first principle studies of graphene/EuO 
hetero-structuresi^. This is due to band folding. 

The lattice constants of the graphene-EuO unit cell 
are exactly three times those of the graphene unit cell. 
Hence, the reciprocal lattice constant of the commen¬ 
surate BZ is ^ that of graphene’s BZ as shown Fig. 
jljb). The outer hexagon (red - online) is the BZ of the 
graphene primitive cell, and the central hexagon (green 
- online) is the BZ of the hetero-structure unit cell. Fig. 
mb) shows that the K and K’ points of the graphene 
BZ lie at equivalent F points in the extended zone of 
the hetero-structure BZ. This results in zone folding of 
graphene’s K and K' points to F. The differences in 
the dispersions and non-linearity are due to inter-valley 
coupling. We address this feature when we construct the 
model Hamiltonian to describe the band dispersions of 
the two hetero-structures in section HI. 


B. Band dispersion with SOC 

Our ab initio calculations that include SOC show very 
different behaviors of the spin resolved bands in the two 
hetero-structures. Since the spin-resolved bands intersect 
in the misaligned hetero-structure, the addition of spin- 
orbit coupling in this system will be more pronounced 
than in the aligned hetero-structure. In the aligned 
hetero-structure spin-orbit coupling leads to a small split¬ 
ting of the spin resolved bands and the dispersion remains 
gapped, therefore, we focus on the effect of spin-orbit 
coupling in the misaligned hetero-structure. 

In Fig. El we plot the band dispersion of the Eu- 
misaligned hetero-structure along path F — K. Fig. [3] 
shows that SOC breaks the degeneracy of the bands 
shown in the left inset of Fig. [5] and gaps the bands. 
The conduction bands and the valence bands split by 8 



FIG. 3: (Color online) Band structure with SOC of 
Eu-misaligned structure calculated along the path F to K 
where K is 0.57 A“^ away from F. 
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Structure 

Eg imeV) 

A.|- {meV) 

A 4 , (meV) 

Se (meV) 

Sh (meV) 

Eu aligned 

127 

309 

344 

182 

217 

Eu misaligned 

-38 

173 

182 

211 

220 


TABLE I: Energy gaps of the EuO-graphene-EuO structures at the Dirac point. Eg is the bandgap of the gapped Dirac 
cone. A-|- is the spin-up gap, and Aj_ is the spin-down gap. The spin-splitting of the electron and hole bands at E are 5^ 
and 5h, respectively. 


meV and 26 meV at the F point, respectively. There are 
two local minimum gaps between the conduction band 
and valence band in the band dispersion near F with val¬ 
ues of 0.1 and 0.3 meV, respectively. The gap between 
the two conduction bands is 0.2 meV; while the gap be¬ 
tween the two valance bands is 1.2 meV. The resultant 
bands show interesting topology on either side of the F 
point which we address using a model Hamiltonian in the 
next section. 


III. LOW-ENERGY EFFECTIVE 
HAMILTONIAN 

In graphene, the gapless Dirac cones at K and K’ 
are protected by time-reversal and inversion symmetry. 
Since these Dirac points are separated in the BZ, small 
perturbations cannot lift this valley degeneracy. There¬ 
fore, the valley index is a good quantum number. In the 
?>N X ZN unit cell, due to zone folding of graphene’s BZ, 
both valleys K and K’ get mapped to F. Hence, valley 
symmetry is no longer preserved and inter-valley interac¬ 
tions can gap the Dirac bands at F without breaking in¬ 
version or time-reversal symmetry. In the graphene/EuO 
hetero-structures, Eu adatoms positioned at the bridge 
and hollow sites contribute two distinct inter-valley in¬ 
teraction terms that are responsible for the non-linear 
dispersions obtained from the ab initio calculations. In 
this section, we construct a model Hamiltonian that cap¬ 
tures the effect of these inter-valley interaction terms, 
and we analyze their effect on the band dispersion. 

A. Inter-valley interactions 

The following model Hamiltonian that acts on an 8 
component spinor is consistent with the lattice symme¬ 
tries, and it describes the salient features of the band 
dispersion near the F point in the absence of SOC. 

iJo = vf{^xTzPx + o-yPy) + AexSz + Ayfx + ma^fx (1) 

In Eq. o, Ti, Gi and Si are the standard Pauli matrices 
acting on the valley, sublattice, and spin degree of free¬ 
dom, respectively. The first term is the standard low- 
energy Hamiltonian describing the linear dispersion of 
the Dirac bands in graphene at the two valleys = ±1 
that are now folded to F. The second term is the ex¬ 
change coupling term induced by the magnetic moment 
of the Eu atom resulting in proximity induced exchange 


splitting Aea: between the spins. The last two terms of 
Eq. ([T} capture the influence of the Eu atoms on the 
graphene layer. In both the hetero-structures of Fig. 
[IJc,d), Eu atoms can sit on a C-C bond, referred as 
the bridge site, and in the middle of the hexagon, re¬ 
ferred as the hollow site. The position of the bridge 
Eu atom reduces the graphene lattice symmetry from 
t Civ resulting in the term A^Tx. in Eq. (P). This 
term corresponds to a valley pseudospin Zeeman term in 
x-direction^ and shifts the Dirac cones from F = (0, 0) 
to (0,±Ai,/rn). The last term, m&xfx results from the 
Eu atom sitting at the hollow site of graphene hexagon; 
we refer to it as an inter-valley scattering term. The 
combined result of these terms, along with the relative 
strengths of A^x, Ay, and m, give a rich band dispersion 
and also account for the differences in the band disper¬ 
sions of the two hetero-structures that we explore next. 

The difference in the band dispersions of the two 
hetero-structures is related to the relative magnitudes of 
Aexi Ay, and m. The energy dispersion of model Hamil¬ 
tonian Hq is 

E± = ±Aex ± _|_ t,|,|p|2 -I- A2 ± 2 A„y^ -I- v^pI, 

( 2 ) 

where |p| = ^Jpx + P^- For Aex = 0 the band dispersion 
has two important features, if m > A„ the dispersion is 
elliptical and gapped at F {px = Py = 0), with an en¬ 
ergy gap 2\Ay — m\. In contrast when A^, > m the Dirac 
points shift from F to (0, i-^/A^ — w? jvp) and graphene 
retains its semi-metallic structure with two Dirac cones 
at ( 0 , A^ — w? jvF)- For Aex ^ 0 and m > A^, there 

are three possibilities determined by the relative magni¬ 
tudes of Aex and |A„ — m|. When (a) Aex < \m — A„|, 
there is a clear gap between the spin resolved states in 
Fig. Pa). For (b) Aex > 1”^ — A„|, the band dispersion 
exhibits an overlap between spin-up and spin-down bands 
shown in Fig. Pb). Finally, at the transition point be¬ 
tween scenarios (a) and (b) when (c) Aex = |ni —A„|, the 
elliptical bands touch. This indicates that when m > Ay 
there a is critical value of Aex > \m — A„| at which 
the spin resolved states intersect. On the other hand, if 
Aex ^ 0 and A^, > m, the shifted Dirac points which 
now appear at (0, ±a/A^ — m?/ vf) exhibit crossing of 
spin-resolved bands indicating that spin-resolved bands 
cross for any value of Aex 0. The gap at F also de¬ 
pends on the relation between Agx and | A„ — m\. When 
Aex < |A„ — m|, the band structure is shown in Fig. 
Pc), and the case of Aex > |A„ — m| is shown in Fig. 
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FIG. 4: (Color online) Band dispersions in the absence of SOC for different values of m, Asx, and A„. (a) m > A„ and 
Aex < |A„ — m|, (b) m > A„ and A^x > |A„ — m|, (c) m < A„ and A^x < |A„ — m|, and (d) m < A„ and A^x > |A„ — m\. 


Relation between m and A„ 

value for Aex 

band dispersion 
along Px 

band dispersion 
along Py 

m > Av 

Aex < m — A„ in Fig. (Ua) 

parabolic, spin 
splitting 

parabolic, spin 
splitting 

m — Av < Aex in Fig. Ulb) 

parabolic with an 
overlap between 
spin-up and 
spin-down bands 

parabolic with an 
overlap between 
spin-up and 
spin-down bands 

m < Av 

Aex < |m — A„| in Fig. (He) 

parabolic with spin 
splitting 

two Dirac cones with 
spin splitting 

Aex > m - A„ in Fig. lUd) 

parabolic with an 
overlap between 
spin-up and 
spin-down bands 

a larger overlap 
between the spin-up 
and spin-down 
double Dirac cones 


TABLE II: Description of band dispersions in the absence of SOC for different parameters of Eq. 0. Plots of the 
dispersions corresponding to different relative strengths of the model parameters are shown in Fig. 3] 


lid). Now that we have established the conditions for 
the intersection of spin resolved bands, we explore the re¬ 
sults of spin-orbit coupling on the graphene/EuO hetero¬ 
structures. 


B. Spin Orbit Coupling 

Spin-orbit coupling introduces two additional terms 
consistent with the lattice symmetries, 

^SOC — dySa;) -t- XlO'^Tz- (3) 


that breaks time reversal symmetry. Therefore A/ 7 ^ 0 
for both structures. Since in-plane inversion symmetry is 
preserved in both structures, we neglect the Dresselhaus 
spin-orbit coupling. Our calculations indicate that the 
strength of the spin-orbit coupling represented by Xr and 
A/ is always smaller than m. Ay and A^x, so we restrict 
our discussions to this case. For Xr, Xj < m. Ay, Aex, 
the spin-orbit coupling introduces gaps between spin- 
resolved bands whenever they intersect (for example see 
Fig. 0] (b),(c) and (d)). With the addition of spin-orbit 
coupling the band dispersion becomes gapped, and the 
bands are a linear combination of spin-up and spin-down 
states. 


The first term is the Rashba spin-orbit coupling which 
breaks inversion symmetry in the plane perpendicular to 
the graphene sheet. Hence, A^ = 0 in the aligned struc¬ 
ture. The second term is the intrinsic spin-orbit term 


By fitting the band dispersion with spin-orbit coupling 
obtained from the ah initio calculations shown in Fig. 
[21 we determine the best ht parameters for our model 
Hamiltonian. This gives hvp = 3.5 eV- A, A^x = 80 
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FIG. 5: Band dispersion with SOC of Eu-misaligned 
structure calculated from model Hamiltonian H along the 
path r to K using parameters obtained from ab initio 
calculations. Hvf = 3.5 eV- A, Aei = 80 meV, m = 48 meV, 
A„ = 17 meV, \r = 5 meV, and A/ = 1 meV. 


(a) 


occupied 


:z 


U 


A/0, m=0 

total 


3 


unoccupied 



0.02 


0.04 

AJeV) 


0.06 


0.02 


0.04 
A, (eV) 


0.06 



meV, TO = 48 meV, A„ = 17 meV, = 5 meV, and 
A/ = 1 meV. The band dispersion along the path F — K 
calculated from the model Hamiltonian H = Hq + i^soc 
is shown in Fig. [S] The model Hamiltonian captures 
all four anti-crossing gaps at about the same position in 
momentum space. Next, we study the topological prop¬ 
erties of these bands obtained from H and calculate the 
Hall conductance for a range of band parameters. 


IV. QUANTIZED ANOMALOUS HALL EFFECT 
IN GRAPHENE/EUO HETERO-STRUCTURES 


FIG. 6: (Color online) The Ghern number calculated as a 
function of Aex for 4 different cases of A„ and m. The red 
open circles show the Chern number of the occupied bands, 
the blue ‘x’ symbols show the Chern number of the 
unoccupied bands, and the black triangles show the Chern 
number of the summation of all bands. The Chern number 
of the occupied bands is 2 for all values satisfying (a) 

A„ yf 0, m = 0 and (b) A„ > m, m yf 0. For condition (c), 
A.„ =0, m yf 0, the Chern number of the occupied bands 
becomes 2 for Aex > ui. For this example, m is chosen to be 
0.04 eV. (d) For A„ yf 0 and m > A„, (in this example 
A„ = 0.01 eV and m = 0.04 eV) the topological transition is 
pushed to a higher value of Aex = 0.12 eV. 


In Ref. d, Qiao et. al. found that ferromagnetic 
graphene in the presence of Rashba spin-orbit coupling 
shows the QAH effect with axy = 2e^/h. First princi¬ 
ple calculations also demonstrated that this QAH phase 
can be engineered by doping with 3d or 5d transition- 
metal atoms or the proximity of a layered antiferroma- 
gentic insulator. In all cases studies thus far, the low 
energy bands are at the K and K' points of the hexago¬ 
nal BZ, and the Hall conductance in the gap is quantized 
(Jxy = 2^ Ih as long as Ay? y^ 0 and A^x 0. In the 
graphene/EuO hetero-structure, as shown in sections H 
and HI, the low-energy bands are no longer at K and K' 
but at r, and inter-valley interactions significantly mod¬ 
ify the band dispersion and hence the topological prop¬ 
erties of the bands. Therefore, we now analyze the effect 
of inter-valley interactions on the topological properties 
of graphene/EuO structures in the presence of in-plane 
inversion symmetry. 

The Hall conductance is calculated from the integral 
of the Berry curvature over the BZ of the occupied bands 
and can be expressed as 

n, ^ Jbz \2.Tr) 

where a corresponds to the band index, Ep denotes the 
Fermi energy, ect(p) is the energy eigenstate, and Qa(p) 


is the Berry curvature of the band. The Berry cur¬ 
vature in terms of the band eigenstates can be expressed 
as 


{ua\dH {^)/dp^\up){up\dH{Tp)/dpj\ua) 

(e;3(p)-e.(p))2 

( 5 ) 

where the Einstein summation convention is used for 
the roman indices i and j, Cij is anti-symmetric ten¬ 
sor, and Ua(p) is the band eigenstate. It is in¬ 
structive to note that in-plane inversion symmetry dic¬ 
tates Qct(p) = Qq,(— p) and time reversal symmetry im¬ 
poses Oa(p) = —^2a(—p). For graphene/EuO hetero¬ 
structures time reversal symmetry is broken due to ex¬ 
change splitting caused by the ferromagnetic substrate, 
however in-plane inversion symmetry is preserved. We 
take advantage of the in-plane inversion symmetry by cal¬ 
culating Berry curvature in the upper half-plane py > 0 
and multiplying by a factor of 2 to account for the lower 
half-plane Py < 0. The Berry curvature is calculated nu¬ 
merically. Our calculations satisfy that the sum of the 
Berry curvatures over all the bands is zero at every p 
point in the BZ, as expected from Eq. 

The model Hamiltonian H only captures the low- 
energy bands near F and may not be valid over the full 
BZ of the graphene/EuO hetero-structure. The Berry 
curvature calculated using falls rapidly away from 


Qa(p) = Im 
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the r point. This allows us to restrict our calculations of 
the Hall conductance to a neighborhood of F. It is well 
known that when the Fermi energy lies in the gap <Jxy is 
quantized and the Hall conductance at zero temperature 
can be expressed as 

(^xy = 'y ^ Con ( 6 ) 

OL 

where the prime indicates summation over the occupied 
bands, and Ca is the Chern number of the band 
that we calculate for different parameters of our model 
Hamiltonian H. We next discuss the Chern numbers at 
Ep = 0 for the occupied and unoccupied bands. 

The calculations for the Chern numbers were per¬ 
formed for 4 different cases with fixed values for Xp, m, 
and A„ (with Xp < m, A„) as a function of ■ The 
four different cases are depicted in Fig. [ 6 l;a)-(d). Figs. 
m a) and (b) show the results when A„ > m. In this case 
the Chern number is quantized and gives a Hall conduc¬ 
tance axy = 2e^//i for any value of Xp ^ 0. However, 
when A„ < m and m 0, there is a topological transi¬ 
tion as a function of Agx and the Chern number changes 
from 0 to 2 as shown in Figs. inKc,d). For A„ = 0, the 
transition occurs when Agx > m as shown in Fig. He). 
For A„ ^ 0 and to > A^,, the transition is pushed to 
a higher value of Agx as shown in Fig. Hd)- For both 
(c) and (d), to = 0.04 eV. In (d), A„ = 0.01 eV, and 
the transition occurs at Agx = 0.12 eV. Unfortunately, 
we have been unable to find an analytical expression for 
the topological transition for A„ < to. From the cal¬ 
culations, we conclude that the Chern number is 2 for 
the case of A„ > m and that the system undergoes a 
topological transition for A„ < to as a function of Agx- 

The results can be summarized in terms of the phase 
diagram shown in Fig. [T] In this calculation, Xp = 1 
meV, and Agx = 70 meV are constant, and the behavior 
of Chern number is calculated as a function of both the 
magnitude of the valley pseudospin Zeeman term A„ and 
the inter-valley scattering term to. As shown in Fig. [71 
when A„ > to, the Chern number is always 2 giving 
a Hall conductance Uxy = 2e^/h. At small values of 
A„, a more complicated situation occurs in the region 
TO > A„ of the phase diagram, however, for sufficiently 
large values of A„, the phase transition occurs at to > A„ 
for a fixed value of Agx- 

V. CONCLUSIONS AND OUTLOOK 

Using insights from first principle calculations and 
lattice symmetries, we constructed a model Hamilto¬ 
nian to describe commensurate graphene/EuO hetero¬ 
structures. In commensurate graphene/EuO structures 
band folding maps the Dirac cones to the F point of 
the hexagonal superlattice BZ of the combined hetero¬ 
structures. Apart from inducing proximity exchange 
splitting to the graphene bands, the Eu atoms also in¬ 
troduce two distinct types of inter-valley interactions, a 
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m(eV) 

FIG. 7: Phase diagram as a fuention of m and A„ for fixed 
values of Xp = 7 meV and Aea, = 70 meV. The curve 
indicates the phase boundary between a Chern number of 2 
on the left side of the curve and 0 on the right side of the 
curve. 


valley pseudospin Zeeman term and an inter-valley scat¬ 
tering term, whose strengths are captured by two model 
parameters A^, and m. The combined effect of exchange 
and inter-valley interactions results in a non-linear dis¬ 
persion at the F point which is captured by the model 
Hamiltonian. The parameters of the model Hamiltonian 
are determined by fitting to the band dispersion obtained 
from the ab initio calculations. 

Using the model Hamiltonian with Rashba SOC, ex¬ 
change, and inter-valley interactions we calculate the 
band dispersion and the topological properties of the 
commensurate graphene/EuO hetero-structures. The 
inter-valley interactions can significantly influence the 
topological properties of the bands for non-zero Rashba 
{Xp yf 0) and exchange splitting [Agx ^ 0). For A„ > to 
with Ai{, Agx ^ 0, the commensurate graphene/EuO 
heterostructure is a Chern insulator with a Hall conduc¬ 
tance Gxy = 2e^/h, whereas for small A„ with to > A„, 
the phase diagram becomes more complicated and one 
needs a large exchange splitting or Rashba SOC to re¬ 
alize the Chern insulating phase. Our calculations indi¬ 
cate that even in the presence of in-plane inversion sym¬ 
metry, inter-valley interactions can significantly influ¬ 
ence the topological properties of graphene/EuO hetero¬ 
structures. 

For a random incommensurate crystallographic stack¬ 
ing of graphene on EuO, the inter-valley coupling will 
be negligible since the Dirac cones will remain at K and 
K'. However, since any incommensurate stacking will 
break the in-plane inversion symmetry of the graphene 
layer, the model Hamiltonian H would acquire an addi¬ 
tional term Ma^. In this case the topological prop¬ 
erties will depend on the relative strength of M and 
Agap < niin(Afl, Aea;). The system will exhibit a QAH 
effect with a Chern number 2 only if M < Agap. 



To observe the QAH effect in graphene/ferromagnet 
hetero-structures, it is important that the disorder in¬ 
duced broadening E of the bands be smaller than the 
topological band gap Ag^p. The critical temperature re¬ 
quired to observe the QAH effect is proportional to the 
mobility gap defined as Agap —E, which must be positive. 
In order to increase the topological gap Agap, it is im¬ 
portant to have a large Rashba SOC Xr, which is small 
^ 7 meV in our calculations. The Rashba spin-orbit cou¬ 
pling can in principle be enhanced by hydrogenation or 
deposition of heavy adatoms on the graphene surface^^. 
Alternatively, the spin orbit coupling can also be en¬ 
hanced by creating hetero-structures consisting of transi¬ 
tion metal dichalcogenides (TMDs) which possess a large 
spin-orbit coupling (e.g. TMD/graphene/EuO)^. Even 
for negative values of the mobility gap (Agap — E < 0), 
the graphene/ferromagnetic structures will exhibit an un¬ 


quantized anomalous Hall effect. However, in this case 
the anomalous hall effect will be additionally influenced 
by disorder induced extrinsic effects^! like side-jump and 
skew scattering mechanisms which are beyond the scope 
of this study. 
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